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PACS 03.67.Mn - Entanglement measures, witnesses, and other characterizations 
PACS 03.65.Ta - Foundations of quantum mechanics; measurement theory 
PACS 03.65.Ud - Entanglement and quantum nonlocality 

Abstract - The common use in literature of orthogonal measurements in obtaining quantum 
discord for two-qubit states is discussed and compared with more general measurements. We prove 
the optimality of orthogonal measurements for rank 2 states. While for rank 3 and 4 mixed states 
they are not optimal, we present strong numerical evidence showing that they give the correct 
quantum discord up to minimal corrections. Based on the connection, through purification with 
an ancilla, between discord and entanglement of formation (EoF), we give a tight upper bound 
for the EoF of a 2 A/" mixed state of rank 2, given by an optimal decomposition of 2 elements. 
We also provide an alternative way to compute the quantum discord for two qubits based on the 
Bloch vectors of the state. 
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Introduction. — The discovery and analysis of the 
quantum speedup within algorithms where entanglement 
is totally absent [1 has raised considerable interest to al- 
ternative measures of quantum correlations beyond entan- 
glement. In contrast to the paradigm of a tensor product 
structure as the root of classicality of correlations, a new 
paradigm based on the ignorance produced by measure- 
ment has crystallized into several measures of quantum 
correlations, of which perhaps the most widely used is the 
quantum discord j2^. In a sense, it captures the fact that 
unless measurements on party B leave unaffected party A 
in a bipartite states qab^ we cannot really speak of such 
state as being purely classically correlated. 

An intense recent research activity is based on quantum 
discord as a quantifier for quantum correlations for two- 
qubit states. Yet it is typically used in a simplified form, 
where only orthogonal measurements are considered (see 
e.g. [IlIllllHS] ) • There is now a raising concern about the 
possibility that more general measurements might modify 
the value of quantum discord, thus weakening the conclu- 
sions of some recent works. The work by Hamieh et al. 
took a first step showing the sufficiency of projective mea- 
surements (i.e. rank 1 POVM's) for states of two qubits. 
This of course does not demonstrate that two orthogo- 
nal projectors are enough. In spite of the great reduction 
of the complexity of the problem [9 , the optimal projec- 



tive measurement of a qubit can have between 2 and 4 
elements [10] (the case of 2 elements corresponds to or- 
thogonal measurements) and the question of how many 
elements the optimal POVM has is still open. In fact in 
the work by Hamieh et al. [9l only a very particular state 
is studied. Hence, the matter remains unsettled and us- 
ing orthogonal measurements seems to be an unnecessary 
restriction. 

In this Letter, we show that orthogonal measurements 
are sufficient to obtain the discord of rank 2 states of two 
qubits, while for rank 3 and 4 they give a pretty tight 
upper bound. Moreover, given the relationship between 
quantum discord and entanglement of formation for pure 
tripartite states [11], we give a formula for discord based 
purely on the eigenvectors and eigenvalues of the original 
state, valid exactly for rank 2. We also show that the 
entanglement of formation of a 2 TV rank 2 state has 
a tight upper bound given by optimal decompositions of 
2 elements. Finally, we give an alternative formula for 
quantum discord of two qubits states of any rank based 
on its Bloch vectors. 

Quantum discord. — Two classically equivalent for- 
mulas for the mutual information in a bipartite state, re- 
lated by Bayes rule are X{A : B) = H{A) + H{B) - 
H{A,B) and J{A : B) = H{A) - H{A\B), where H{.) 
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is the Shannon entropy and H{A\B) is the conditional 
Shannon entropy of A given B. Their quantum counter- 
parts, however, differ substantiahy [2], the former being 
known as the quantum mutual information: 



I{g) = SigA) + SigB)-Sig), 



(1) 



where S{.) is the von Neumann entropy and Qa^b) ^ire the 
reduced states after tracing out party B{A). It is pre- 
cisely in J7(A : B) where measurements come into scene, 
since the conditional entropy of A given B is the one given 
by measurement outcomes on party B. Though a general 
measurement has to be associated to a generic POVM, 
for the sake of simplicity the community has used the re- 
stricted set of perfect (von Neumann, or orthogonal) mea- 
surements, i.e. : 

J{g){nf)=S{gA)-S{A\{Iif}) (2) 

with the conditional entropy defined as 6'(A|{nj^}) = 
Y.iPiS{QA\nf)i Pi = TrAB(nf ^) and where Qj^\^b = 
lif glif /pi is the density matrix after a complete projec- 
tive measurement ({11^}) has been performed on B. No- 
tice that n^ are orthogonal projectors. Quantum discord 
is thus defined as the difference 



^A:B 



(q) 



min \S{gB) 



S{g)^S{A\{Uf})] 



(3) 



minimized over all possible orthogonal measurements. 
However, more general measurements should be used to 
exhaust the minimization problem, as already stated in 
the seminal papers [2j[3]. Hence the following generaliza- 
tion is required 

QA\uf -^ QA\Ef = T^^B{EfgAB/Pj) (4) 

where the elements of the POVM Ef fulfill J]^- Ef = 1b 



^A:b{q) = min [S{gB) 



S{g)^S{A\{Ef})]. (5) 



Based on the convexity properties of the conditional en- 
tropy S{A\{EP}) Hamieh et al. [9] (see also [H]) show that 
POVM's which optimize discord are extremal or indecom- 
posable, i.e. they cannot be obtained by mixing other 
POVM's; further, it has been shown ^lOj that extremal 
POVM's for qubits are of rank 1 and can have between 2 
and 4 elements. 

Unified picture through purification. — In [11] it 

was shown that given a mixed state gAB and its purifica- 
tion It/jABc) through an ancilla qudit C, the following rela- 
tion between the conditional entropy when B is measured, 
and the entanglement of formation Ep of the subsystem 
AC holds: 



min S{A\{Ef})=EF{gAc) 



(6) 



between the conditional entropy when B is measured, 
and the entanglement of formation Ef of the subsystem 
AC. The minimization of the AB conditional entropy 
over POVM measurements on B is thus equivalent to 
minimization of Ej? in AC over all ensemble decomposi- 
tions. Hence the number of elements giving the optimal 
ensemble decomposition of EF{gAc) coincides with the 
number of elements of the POVM which minimizes 

s{A\{Ef}) m- 

Theorem 1. Given a bipartite mixed state of two 
qubits pAB of rank 2, the optimal measurement giving 
the quantum discord is a 2 element POVM. The elements 
of such POVM are orthogonal projectors. 

Proof. Consider a rank 2 state of two qubits 
with spectral decomposition pab = XliJi <^i|'0i)('0i|- 
Its purification by an ancilla qubit C has the form 
I^ABc) = ELi V^IV^^)K)c5 where \i)c forms an 
orthonormal basis in the Hilbert space of the an- 
cilla qubit. We can also Schmidt decompose this 
state as \^ abc) = YaLi V?i\^) B\(t>i) AC , where 
m = iniii{dB^dAdc) = dB = "^ with dB is the di- 
mension of the Hilbert space of party B and so forth. So 
the partition AC has the form pAC = Zlill ft|^i)(^i| 
and hence is of rank 2. Wootters [13] showed that the 
entanglement of formation of this two-qubit mixed state 
is obtained from an optimal decomposition made up of 
as many elements as its rank^ which in this case is 2, 
which in turn means that the POVM in B that realizes 
such decomposition has 2 elements. Being optimal 
POVM's of rank 1, the 2 elements of such POVM are 
necessarily orthogonal [10]. This can be seen by noticing 
that a rank 1 POVM of 2 elements Ei = ai|(/)i)(0i| 
and E2 = Q^2|02)(</>2| has to fulfill positivity Ei > 
and normalization Ei -\- E2 = 1, which necessarily lead 
to orthogonality of its elements (this is easy to show 
when the elements Ei are written in Bloch form [10] . QED 

Corollary. The quantum discord of a rank 2 state of 
two qubits is given by 

SaMqab) = S{gB) - S{gAB) + S{C{gAc)). (7) 
with 



Qac = tr^ 




^/X^\^^){^J\^\i)c{j\ 



(8) 



where {Aj, |^i)} is the spectral decomposition of gAB, and 
\ic) is any orthonormal basis in He- The function £ is 
given by 

,1 + yr^c^. 



£{C) = h{- 



-), 



where 



h{x) = -xlogax- (1 -a;)log2(l-a;), 



(9) 



(10) 
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and where C{p) is the concurrence of p [T3l, given by 
max(0,/i — h — h — h)^ with li the eigenvalues of the 
hermitian matrix R{qac)i where R{p) = \J ^fpp^fp and 

P= [(Jy®(7y)p''[(7y®Gy). 

Hence, to obtain the discord, diagonalize qab and use 
its eigenvalues and eigenvectors to construct qbc- Obtain 
the concurrence of qbc^ substitute it in £{C) and obtain 
the optimal conditional entropy between A and B. This 
result was implicit in ref. Jlj. 

Quantum discord of rank 3 and 4 states. — Or- 
thogonal measurements do not give the optimal discord for 
higher rank states, as ffist found in fTS" through a coun- 
terexample based on maximally discordant mixed states 
(MDMS) [16]. Their study was limited to 3 element 
POVM's leading to a deviation of ~ 2 x 10~^ with respect 
to orthogonal projectors. However we will give evidence 
that the set of states where 3 and 4-element POVM's are 
needed is indeed small and the improvement in discord 
using Eq. ([5]) is tiny. 

In the case of states qab with rank higher than 2, a pu- 
rification would yield a qubit-qudit system in AC, whose 
optimal decomposition (for Ep) is not known. Therefore 
no analytical tool can help us discriminate how many ele- 
ments build the optimal POVM for the quantum discord 
of states with rank higher than 2. However, it is known 
that for qubits the optimal measurements are given by 
rank 1 measurements with 2, 3 and 4 outcomes. This 
knowledge is based on two facts: i) the conditional en- 
tropy S{^A\{Ef\) is a concave function over the convex 
set of POVM's 0112], hence only extremal POVM's wih 
minimize it, ii) extremal POVM's of qubits are rank 1 
and have 2, 3 or 4 elements [10:, i.e. E^ = aj\(j)j){(j)j\ 
(j = 2, ...A/", A/" < 4), with a real and nonzero and \(j)j) are 
pure states (nonorthogonal unless A" = 2). 

Even considering at most 4 elements POVM's, the nu- 
merical analysis of this problem is challenging, as detailed 
in the following. We start parametrizing the N elements 
of a POVM as: 

El = ai|0)(0| 



E, 



a 



.u{0j,<pjmmHoj,<f>j) {i<j<N), 



where U{6^ (j)) is a qubit rotation, plus a final global 
rotation U{Q.^ ^) acting on all elements; i.e. Ei = 
U{Q,^)EiU\Q,^). The completeness relation Xli ^i = 
1 solves the coefficients in terms of the angles: a^ = 
ai{{Oj^(j)j}). This means running 6(8) loops in angleq^, 
for 3(4)-element POVM's, for each state, in order to solve 
the minimization problem in the discord definition (J5j) . We 
must note that orthogonal measurements are a limit case 
in the definition of 3-element POVM's; in the same way 4 
element POVM's do not include the case of 3 elements. 



The numerical evaluation of discord is very sensitive to 
identification of the minimizing POVM and a proper scan 
of all possible POVM's requires small step sizes for the an- 
gles {Oj,(j)j}. To give an idea on the sensitivity of the mini- 
mization on the step size in the angles, we show in fig.[T]the 
values of discord for 2, 3 and 4-element POVM's ((^2, ^3, ^4) 
against the step size for three states. In all cases it is sur- 
prising the oscillatory nature of discord even for rather 
refined samplings (small angles). This demonstrates the 
importance to scan all POVM's over different step sizes 
to gather the minimum, best approximating Eq. (J5]). We 
notice that a refinement until angular step sizes ~ 0.027r 
is feasible only for 2 and 3 elements POVM's. In any 
case, a good level of accuracy is obtained if the lowest 
value of discord obtained among different angular preci- 
sions is retained. The insets of fig. [1] show the minimum 
value of discord obtained inside a box of angular precisions 
Ai9 = A0G [w,0.257r]. 

In fig. [1^ we show the state with highest deviation we 
have found (highest point in fig. [2]) in a scan of 10^ random 



states of rank 3 and 4, with 62 — S[ 



'3(4) 



10- 



This devi- 



-•^ Though for the case of 3-element POVM's the number of loops 
can be reduced to 5 by using geometrical arguments. 



ation is high above the typical deviation we have found of 
around 10~^. In fig.[D3 we show a MDMS of rank 3, sepa- 
rable, but with maximum discord versus classical correla- 
tions [H]: ^MDMs = (1 - 6)(m|00)(00| + (1 - m)|ll)(ll|) + 
e|^-)(^-| with 1^-) = (|01) - |10))/a/2 the usual Beh 
state and values e = 0.2349602, m = 0.11. In this case, 
POVM's with more than 2 elements are needed. This is 
actually a rather singular event, as we will see in fig.[2l In- 
deed in fig.IIb we show that an improvement of ~ 8 x 10~^ 
is provided by 3, 4-element POVM's (^3(4)), as compared to 
orthogonal measurements (^2)- The most common situa- 
tion is represented in figure [it, for a generic state obtained 
by a random density matrix of rank 3. 

Next we aim to establish the abundance of states for 
which 3 or 4 elements POVM's provide the improvement in 
discord found in fig.[Tl We then start with a scan of Hilbert 
space where 10^ random density matrices of rank 3 and 4 
have been generated according to the Haar measure. We 
plot in fig. [2] the deviations 62 — ^3(4) (only when positive) 
versus the result for orthogonal measurements (62). In this 
figure we have sampled the angles from steps ~ 0.37r until a 
lower limit of w= 0.037r (O.IStt) for 3 (4)-elements POVM's 
respectively. We observe that the optimal discord is given 
by orthogonal POVM's, except for a 0.63%(0.001%) of 
states which have a typical deviation of order 10 ~^ (see 
Table [T]) . In order to discriminate the dependence of such 
abundance and degree of deviations on the scan size (A^), 
or the precision in the angles (A^ = Acj) G [w,0.257r]), we 
present in Table [T] the results obtained for different sample 
characteristics A^,w. 

We find that for a given angular precision (for instance 
up to w= 0.037r in Table [T]) the abundance of deviant 
states (p) is not sensitive to the sample size (the cases 
A^ = 3 X 10^ and 10^ can be compared in the table). On 
the other hand discriminating states with 62 > Ss in which 
these values are close, requires a quite refined angular pre- 
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Fig. 1: (color online) Quantum discord minimized by 2 (^2, 
black), 3(^3, red) and 4(^4, orange) elements POVM's with the 
given step size A^ = A(/) in the angles of the POVM elements 
Ei. The states are a) the most deviant state we have found 
(highest point in fig. 0, b) a MDMS of rank 3 with m = 0.11, 
see [16], and c) a random rank 3 state. In the inset we show 
the precision of fc, ^3 and ^4, by plotting its minimum inside 
a box A6' = A0 e [w, 0.257r]. 



cision (w), therefore leading to a highly dependent abun- 
dance (p) (see strong variation of p to find states with 
62 > Ss for N = 10^ in the table). Due to computational 
limitations such a deep numerical study was only possible 
for 3-el. POVM's while for 4-el. POVM's we were only 
able to reach precisions up to w= O.IStt, for which we 
have found only one deviating state. At the best precision 
reached we see that the typical (average) deviation 62 — Ss 
is of the order 10~^ with a standard deviation of order 

The set of states for which we find improvements with 
^3(4) is rather small. An interesting question is whether 
these states lie in the neighborhood of rank 3 MDMS [151 
[T6] or are distributed everywhere in the {J'^X} diagram 
(figure 1 in the same reference). We first stress that the 
state in fig. [1^ is nowhere near the MDMS border (it has 
{J c^ 0.17, X c^ 0.1809}). Moreover we can gain some 
insight about the smallness of the MDMS neighborhood 
by investigating the state of rank 3 in fig. [T]d [16] whereby 
we perturbate it with a Bell state. That is we study the 
state ^ = (1 - A)^MDMs + A|<l>+)(<l>+| with |<^+) = (|00) + 
|ll))/>/2. The effect of this perturbation is to move the 
state away from the border in the {J7,X} diagram. We find 
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Fig. 2: (color online) Deviation of quantum discord fe — ^3(4) as 
given by 3- (blue) and 4-element (red) POVM's, for a scan of 
10^ states of ranks 3 and 4. The density matrices are randomly 
generated following the Haar measure. The only state found 
with S2 — S4 > corresponds to that of fig. [T^, with highest 
deviation. 



that already for A c^ 0.002 we reach the transition where 
^2 = ^3 = ^4, meaning that, at least in the neighborhood 
of the MDMS border, 3,4-element POVM's are needed in 
a very tiny region. 

Summarizing our numerical analysis for randomly gen- 
erated states, we have found that orthogonal measure- 
ments are almost enough in the following sense. The im- 
provements ^3(4) < 62 do occur but they represent small 
corrections (a maximum deviation 10~^ was found for only 
one state in a sample of 10^ while other deviations were up 
to 10~^, fig. [2]) and they can be numerically appreciated 
only for a tiny set of states (with the mentioned mini- 
mization up to precision w= 0.027r these states appear 
with probability ^ 10~^ for 63, while for S4 they appear 
with probability ~ 10~^ at w= O.IStt). 

Given this numerical evidence we give the following up- 
per bound: 

Observation. — The entanglement of formation of a 
2^ N rank 2 bipartite state qbc has a tight upper bound 
given by an optimal decomposition of two elements: 



Ef{qbc) < min ^pfc^(|(/)|^)), (11) 

with E{\^%c)) = S{tTB{\cl)%c){^Bc\))' The deviation 
from the equality is on average of order 10~^. 



To show this, take the 2 > 



form Qbc = ELi ^' '''^'^ 



>^k\^Bc){^BC 



N^ rank 2, state in spectral 
It can be purified by 



an ancilla qubit as IV^abc) = J2k^i 
\e\) any orthonormal basis in I-La- 



/^I4>l^lc)' with 
Again according to 



EpigBc) = min Si{Ef}\B) < min S{{Ilf}\B), (12) 
{Ef} {nf} 
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Table 1: Probability (p) to find states with A = ^2 — ^3(4) > in 
a sample of N states whose discord is calculated with angular 
precision up to w. We show the average deviation (A) and its 
standard deviation cr(A). The standard deviation is reduced 
by an order of magnitude if we remove the state a) of figlU We 
do not show the average deviation for 62 > S4 nor its standard 
deviation since only one state was found in the full scan. 
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where restriction to orthogonal measurements 'spoils' the 
minimization. Given the numerical evidence we have pro- 
vided, the improvement of doing full minimization is on 
average at the level of 10~^ (see Table [T]). 

The conditional entropy for a given orthonormal mea- 
surement n^ = |^^)(^^| is given by 

2 2 

S{{n^}\B) = Y.PkS{p%) = Y.PkEip'hc) (13) 

with pk = tr(n^^ABc) and p% = tTAc0^toABc)/Pk' In 
the last equality we have used the fact that p^^ is pure. 
This can be seen by writing explicitly 

P%c = trA(n^^ABc) = 



\^i3c){^%c\ 



'bc\ 



(14) 



with \(j)i^c) = Ei=iei(fc)\/Al|V'Bc). and ei{k) = {^''\ei). 
So, finally, S{p%) = E{\4>'hc))- QED 

We stress the fact that instead of minimization over en- 
semble decompositions with a number of elements ranging 
from R to R^ {R is the rank of the state), as shown to be 
sufficient by Uhlmann [T7], we can safely restrict to de- 
compositions with 2 elements if we are not interested in 
states which are rare to find and have deviations which are 
probably quite small, as we have seen when perturbing a 
MDMS extremal state. 

Bloch formula for quantum discord. — We finish 
by giving an alternative formula for the computation of 
quantum discord of generic two-qubit states of any rank. 



Writing the POVM elements in Bloch form 

Ef =ai{lA^ni -(Ja), 



(15) 



with the positivity and normalization (completeness) con- 
ditions 

<^i > , ^Qfi = 1 , ^aifii = , {\fii\ = 1), (16) 

i i 

and the density matrix qab (coming from purification of 
qbc plus tracing subsystem C) given also in Bloch form 

QAB = - itAB^a-aA^b-as^'Ycia'A^a^ I (17) 

(notice that when we write E-ga^ we mean a- a a® 1b), 
we enunciate the following theorem: 

Proposition. The discord 5a:b{qab) with qab of any 
rank^ written as in eq. (fT7|) . is given by eqs. (|18ti2"T]) and 
m{= 2,3,4) is the number of elements of the extremal 
POVM. 



The minimization is restricted by the conditions for ex- 
tremality [10] for each number of POVM elements (m): 



m=2: 



m=3: 



m=4: 



all extremal, from normalization they follow: 

a^ = - , ni = — n2 = n (i.e. orthogonal) (22) 



all extremal, from normalization they follow: 
as = 1 - ai - a2 , fis ■ 



— (ni+ns) (23) 
as 



fii not in the same plane; 

4 
normalization yields >, <^i^i 



(24) 



Proof. Simple algebra yields the probabilities and out- 
comes of each measurement: 



Pk 



ttB {EfQAB) IPk = 



lilA 



{a + chk) ■ a A 

1 + 5 - Hfe 



(25) 



(26) 



whose entropy can be calculated from the eigenvalues of 
p\ {X^{a^b^cnk)) as defined in eq. (p!9|) . The condition 
that for m = 4 extremality implies that the POVM 
elements cannot lie in the same plane, was derived in [10], 
as well as the normalization conditions. 

In summary, we have proven that the quantum discord 
of rank 2 two-qubit mixed states is obtained using only 
orthogonal projectors as measurements. Strong numerical 
evidence has been given to conjecture that they are almost 
sufficient for higher ranks, except for states which appear 
with probability '-«^ 10~^ and have negligible deviations on 
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Sa:b{QAb) 

\f(a,h,chi) 

Sb 
Sab 



Sb-Sab^ min y^ ai{l ^ Hi -b)y^ H{Xf{a,b,cni)) 



1± 



a + chi 



l^b-fii 
S{gAB) 



C^i i^^x^i^X") ^y^i,yi ^z^i,z j 



(18) 

(19) 

(20) 
(21) 



average of order 10~^ in a sample of 10^ states (w= 0.027r 
for ^3). We discussed the importance of the states sample 
size and of the precision in scanning all possible measure- 
ments, showing the need of using a minimization proce- 
dure over different angular step sizes. Two examples were 
given in figs. [1^ and [TJd, both of rank 3, showing that 3 
and 4 element POVM's give a better quantum discord, 
though with a very small improvement (a maximum of 
order 10~^). Based on the connection [11] between con- 
ditional entropy and entanglement of formation we have 
given a related tight upper bound, namely that the en- 
tanglement of formation of a 2 (g) A^ system of rank 2 
is obtained by 2 element decompositions with very high 
probability and precision. Finally, an alternative quan- 
tum discord formula for generic mixed states of two qubits 
was given in terms of the Bloch vectors of the state, where 
minimization is performed over the Bloch form of POVM's 
elements. As a side remark, we note that Gaussian discord 
in continuous variable systems is also an example where 
projective, but not orthogonal measurements are optimal 

[Sill]. 
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